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Section  1 

PLASMA  THEORY 


A.  F I LAMENTATION  IN  CHARGED  PARTICLE  BEAMS 
Lee  Buchanan 

Since  the  classic  paper  of  Weibel  (1959),  It  has  been  known 
that  homogeneous  plasmas  with  anisotropic  velocity  distributions  are 
unstable  to  transverse  electromagnetic  modes.  In  a streaming  plasma, 
these  modes,  also  known  as  "tearing  modes",  manifest  themselves  as  the 
transverse  pinching  of  the  beam  into  small  current  filaments.  This  fl la- 
mentation results  in  a significant  transverse  heating  of  the  beam,  inducing 
radial  expansion  which  grossly  affects  propagation,  focusing,  etc.  A 
simple  example  is  given  by  Benford  (1973)  who  considers  an  infinite,  homo- 
geneous beam  streaming  with  a velocity  vq  through  a collisionless  plasma. 
The  beam  is  taken  to  be  charge  neutralized  by  the  plasma  and  the  transverse 
velocity  distribution  is  taken  to  be  Lorentzian  (allowing  closed  form 
solutions).  The  dispersion  relation  is  found  to  be 
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where  k is  the  transverse  wave  number.  When  v,  < v there  are  unstable 
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modes  with  maximum  growth  rate 
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clearly  showing  the  stabilizing  effect  of  vq.  This  can  be  interpreted  as 
the  tendency  for  an  anisotropic  beam  to  seek  isotropy  using  the  mechanism 
of  the  t nstabi I i ty . 

Exhaustive  treatment  of  both  the  linear  and  non-linear  devel- 
opment of  the  instability  is  found  in  the  literature.  Rather  complete  non- 
linear analytic  descriptions  of  the  Weibel  instability  are  compared  to  one- 
and  two-dimensional  simulations  by  R.  C.  Davidson  et  al.  (1972)  and  R.  L. 
Morse  and  C.  W.  Nielson  (1971).  Of  particular  interest  are  the  two-dimen- 
sional simulations  of  R.  Lee  and  M.  Lampe  (1973)  of  the  f i I amen tat  ion  of 
relativistic  electron  beams  in  which  the  microstructure  is  described. 

In  the  regime  where  the  background  plasma  is  highly  colli- 
sional  and  conductivity  is  low,  the  f i I amen tat  ion  is  non-res i st i ve  in  char- 
acter and  is  more  generally  termed  fluting  (interchange)  and  firehose.  The 
two  have  received  quite  adequate  consideration. 

In  the  case  of  high,  but  finite,  conductivity,  E.  P.  Lee  and 
F.  W.  Chambers  have  found,  for  a fully  charge  and  current  neutralized 
beam  in  a highly  ionized  background  of  conductivity  o,  the  dispersion  relation 
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where  k^  and  k|  are  the  perpendicular  and  parallel  wave  numbers  and  vq  is 
the  beam  velocity  as  before.  Substitution  of  p = enu  where  u-e/m(v-  iw) 
and  letting  the  collision  frequency  v and  the  parallel  k go  to  zero  recovers 
the  Benford  result.  For  finite  p.  Chambers  (1978)  finds  that  for  a self- 
pinched  beam  (e.g.  Bennett,  1939)  one  can  invoke  the  relation 
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and  the  condition  that 
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where  t is  the  beam  pulse  duration  and  t is  the  magnetic  decay  time  (t 
p m 3 m 
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■ Anp/c  k^),  to  obtain  the  total  growth  for  a transverse  perturbation 
to  a beam  streaming  in  the  z-direction.  This  can  be  written  in  terms  of 
a parameter 
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so  that  one  obtains  as 


total  growth 
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This  analytic  result  is  based  on  the  assumption  of  constant,  isotropic  con- 
ductivity of  the  background  plasma.  Since  the  growth  is  a sensitive  compe- 
tition between  a growth  term  and  a damping  term,  both  functions  of  p,  one 
may  conclude  that  the  t ransverse  structure  of  p and  particularly  the  time 
evolution  of  p due  to  the  Ionization  of  the  background  by  the  beam  may 
be  crucial  to  the  onset  and  development  of  the  instability. 

Presently  we  are  undertaking  to  simulate  in  one  and  two 
dimensions  the  Weibel  instability  in  charged  particle  beams  in  a background 


of  finite  conductivity.  We  intend  to  construct  a particle/fluid  hybrid 
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simulation  code  in  which  the  beam  consists  of  finite  sized  particles  and 


the  background  is  taken  to  be  a fluid  plasma  of  finite  conductivity  o. 
It  is  expected  that  the  new  physics  should  appear  as  a consequence  of 
both  the  self-consistent  time  and  spatial  development  of  the  fluid  con 
ductivi ty. 
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A. 


DRIFT-CYCLOTRON  INSTABILITY  P ARTICLE  SIMULATIONS 
Jae  Koo  Lee  (Prof.  C.  K.  Birdsall) 


J.  INTRODUCTION 

The  drift-cyclotron  instability  is  a collisionless  instability 
occurring  in  any  finite  plasma  with  a Maxwellian  velocity  distribution  in 
a magnetic  field.  This  instability  might  constitute  a part  of  the  major 
plasma  confinement  problems  in  e.g.,  most  mirror  machines  like  2 X II B or  TMX 
even  after  the  dr i ft-cyc lotron- loss-cone  mode  is  stabilized. 

Recent  theoretical  studies  jl)  of  the  dr i ft -eye lotron  instab- 
ility have  attracted  much  attention;  however,  they  have  had  some  diffi- 
culty in  providing  useful  information  on  the  nonlinear  evolution  of  this 
instability.  Our  object  is  to  use  particle  simulation  to  obtain  linear 
and  nonlinear  behavior  of  this  instability. 

One  recent  particle  simulation  of  the  drift-cyclotron  instab- 
ility is  one- dimensional  [2],  treating  the  electron  effects  as  those  of 
a linear  fluid  and  assuming  a fixed  density  gradient,  thus  reducing  the 
spatial  dimension  in  simulation  from  two  to  one.  These  results  provide 
useful  information  but  have  the  shortcoming  of  a fixed  density  gradient. 

We  have  simulated  this  instability  using  a two-dimensional 
particle  code  fer  a Maxwellian  velocity  distribution  with  a density  gradient 
in  a uniform  magnetic  field.  During  the  linear  stage,  our  simulations  show 
exponential  growth  over  a few  ion  cyclotron  periods.  We  compare  this  linear 
behavior  with  linear  theory  which  does  not  use  the  local  approximation 
because  a./R  is  large  (here,  a.  Is  the  ion  gyroradius  and  R is  the 

ip  i p 

density  scale  length).  This  nonlocal  theory  is  that  of  ROOTS  [3|.  which 


has  a periodic  guiding  center  density  profile  unlike  the  bounded  one  in 
our  simulations.  The  linear  growth  rates  and  (in  less  degree)  the  fre- 
quencies in  the  simulations  are  found  to  agree  well  with  the  nonlocal 
theory  predictions  (and  poorly  with  the  local  theory).  The  comparison  of 
simulation  with  the  nonlocal  ROOTS  predictions  is  given  in  Sec.  Ill  for 
two  mass  ratios,  i.e.,  mj/m^^lOO  and  25. 

At  saturation,  the  total  electrostatic  field  energy  reaches 
a few  percent  of  the  initial  ion  kinetic  energy  for  most  runs.  These 
simulation  saturation  levels  are  compared  with  some  existing  nonlinear 
theories  in  Sec.  IV.  At  and  after  saturation,  the  phase  space  plots  of 
electrons  and  ions  show  vortex-like  structure  corresponding  to  the  most 
unstable  mode;  the  density  profiles  of  each  species  become  flatter  as  has 
been  observed  in  lower-hybrid-drift  ins.ability  simulations  [4). 

In  Sec.  II,  a brief  description  of  our  simulation  model  is 
given.  Growth  rates  and  frequencies  are  compared  with  theoretical  pre- 
dictions In  Sec.  Ill  for  two  cases.  In  Sec.  IV,  nonlinear  aspects  of  this 
instability  are  discussed  Including  comparison  with  nonlinear  frequency 
shift  theory,  trapping  theory,  density-gradient  modification,  and  electron 
heating.  Section  V is  the  conclusion. 

II.  SIMULATION  NOVEL 

For  our  study  of  the  density-gradient-driven  drift-cyclotron 

instability,  we  need  to  use  a fully  two-dimensional  particle  code  for  both 

electrons  and  Ions,  I.e.,  the  evolution  of  four  phase  space  variables  in 

the  plane  perpendicular  to  the  magnetic  field  (x,y,v  ,v  ) of  all  the  parti- 

x y 

cles  needs  to  be  followed.  The  common  guiding  center  approximation  for 
electrons  was  not  used  in  our  simulations;  electrons  and  ions  are  both 
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magnetized  and  fully  nonlinear.  The  code  we  use  is  a two-dimensional  electro- 
static one,  EZOHAR  |5).  The  coordinates  and  the  density  profile  of  the 
plasma  slab  model  are  shown  in  Fig.  I.  For  simplicity,  all  the  simulation 
runs  in  this  paper  have  no  temperature  gradient  and  Initially  cold  electrons. 

The  initial  particle  spatial  loading  is  as  follows:  first, 
the  ion  guiding  center  positions  are  determined  from  the  following  profile, 


n.  (x  ) - 
' 9C 


1 + cos  < 

a 


x < n a , Ac 
gc  i 


x > na,/ic 
gc  - I 


Here,  x ■ x + v /u>  . , k i s a constant  indicating  the  density  scale  lenqth 
gc  y c i * 

(i.e.,  <c»a./R^),  and  is  the  ion  cyclotron  frequency.  Then  the  ion  gyro- 
radius  is  added  to  these  guiding  center  pos i t ions  to  determi ne  the  actua I ion 
positions  (velocities  are  chosen  from  a Maxwellian  distribution).  If 
this  process  is  repeated  for  electrons,  then  appreciable  charge  separation 
results,  because  the  guiding  center  profile  is  different  from  the  particle 
density  profile  for  the  species  with  large  gyroradius;  these  two  profiles 
are  the  same  for  cold  electrons  as  noted  in  Naitou  et  al.  [6|.  Instead, 
the  actual  positions  of  ions  are  copied  in  our  simulations  to  give  the 
actual  electron  positions,  thus  resulting  in  no  charge  separation  at  the 
initial  loading.  (As  the  instability  developed,  we  sometimes  observed  char^,? 
separation  producing  appreciable  electric  field  along  the  density  gradient 
direction  responsible  for  some  collisionless  plasma  transport.  This  is  dis- 
cussed in  Sec.  IV  (3).]  The  velocity  loading  is  by  "ring  and  spokes"  in  v^ 

Vy  space  with  the  ring  radius  chosen  randomly  to  fit  a Gaussian  and  the 
spokes  (angles)  are  chosen  randomly,  uniformly. 


Along  x,  the  boundary  conditions  in  EZOHAR  are  those  of  Naitou 
et  al.  [6]  at  the  low  density  side  and  those  of  inversion  symmetry  of 
Nevins  et  al.  [7]  at  the  high  density  end,  x«0.  Along  y,  the  system  is 
periodic  due  to  the  slab  model  (cf.  Fig.  1). 

Typical  parameters  are  as  follows  (these  are  to  be  varied 
widely  but  used  unless  otherwise  specified  in  later  Sections) 

number  of  spatial  grids  * 32  * 32 

number  of  electrons  or  ions  * 16  * 16  x 32  **  8192 

uj  At  'V  0.5 

pe 

k s (a.dn/dx/n)  'v  0.5 

i x *Lx 

m./m  - 25-300 

L « L - 8a. 

x y i 

A.,  -v  0.2a. 

Di  i 


where  <*>  is  the  electron  plasma  frequency,  At  is  the  time  step,  L and  L 
p®  x y 

are  the  simulation  dimensions  in  x and  y,  and  A^.  is  the  ion  Debye  length. 


III.  CHECK  WITH  LINEAR  THEORY 


(1)  m./me«25  Case 


In  order  to  check  the  growth  rates  and  the  frequencies  pre- 
dicted by  the  linear  nonlocal  theory  given  by  ROOTS  [3],  we  made  an  all- 
mode  run  (i.e.,  no  artificial  truncation  of  Fourier  modes  was  used).  The 
parameters  are  (in  addition  to  those  in  Sec.  II) 


m./me  * 25 


k - 0.5 

b>  /u)  * 1 

pe  ce 


2 2 
thus  in  ./it)  . 
pi  c i 
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The  squared  potentials  for  modes  5 and  9 are  shown  in  Fig.  2A  and  2B,  which 
tend  to  follow  the  theoretical  growth  rates  drawn  as  straight  lines.  Reason- 
able agreement  in  growth  rates  extends  to  most  graving  modes  in  the  range, 

1 . 2 < ka . < 3 • 6 as  summar i zed  in  Fig.  3B , which  does  not  i nc lude  the  finite 
grid  (Ax)  effect.  We  observed  less  satisfactory  agreement  in  the  measure- 
ment of  oscillation  frequencies.  The  solid  lines  in  Fig.  3A  show  the  non- 
local ROOTS  theory  [ 3 1 - The  simulation  frequencies  are  measured  during 
approximately  the  first  two  ion  cyclotron  periods  (i.e.,  0 < . t < 1 2 as 

in  Fig.  3)  which  is  considered  to  be  the  linear  stage.  Figure  A is  a fre- 
quency spectrum  for  mode  3 (i.e.,  ka.*1.2  in  Fig.  3)  showing  a quite  broad 
spectrum  with  the  width  of  about  J-l  of  the  peak  frequency.  This  broad 
width  is  due  to  the  short  time  span  used  in  fast  Fourier  transform;  it 
could  be  reduced  by  having  a longer  linear  growth  period  with  just  the 
growing  mode  excited.  The  frequencies  of  modes  7 and  8 in  Fig.  3A  also 
had  a broad  width  around  the  peak  u>“0,  indicating  the  possibility  of 
a vortex  mode. 

(2)  m./m  * 100  Case 
1 e 

Another  check  with  the  linear  theory  was  sought  in  a simulation 

run,  where  larger  mass  ratio  (i.e.,  m./me  = 100)  and  larger  density  gradient 

2 2 

(i.e.,  k * 0.  7)  were  used  in  a lower  density  plasma  (i.e.,  ^pe^ce  ~ 0. 15)  • 

The  growth  rates  are  shown  in  Fig.  5;  agreement  is  again  fairly  good  consi- 
dering that  all  the  growing  modes  in  this  range,  1 < ka . < 6,  wera  al  lowed  to 
grow  in  a single  run.  The  real  frequencies  produced  broad  spectra  as  in  the 
previous  small  mass  ratio  case. 
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II/.  NONLINEAR  BEHAVIOR 


For  the  study  of  the  nonlinear  evolution  of  the  d. ; t-cyclotron 
instability,  our  simulation  uses  single-mode  runs;  namely,  only  one  mode  in 
the  y-direction  is  allowed  to  vary  while  the  perturbed  quantities  (charge 
densities  and  potentials)  corresponding  to  all  the  other  mode  numbers  are 
truncated  at  all  times  after  Fourier  transform  in  the  y-direction.  This 
truncation  includes  the  mode  k^ * 0 which  implies  charge  neutrality; 
independent  of  this  truncation,  s i mu lat ions  star  t exact ly  charge-neutral. 

This  single-mode  assumption  has  been  also  used  in  the  nonlinear  theory  of 
a beam-plasma  instability  by  O'Neil  et  al.  [8],  and  is  a reasonable  one 
especially  in  an  instability  study  with  a discrete  mode  spectrum. 

Since  exact  particle  movers  are  used  for  electrons  as  well  as 
for  ions,  our  simulations  are  restricted  to  using  small  mass  ratios  (e.g., 
m./m^  * 25-300)  to  keep  the  simulation  errors  and  cost  low  when  reasonably 
large  time  increments  are  used.  This  error  or  numerical  heating  was  1-2% 
of  the  initial  ion  kinetic  energy  up  to  the  saturation  time  t,  which  is 
defined  as  the  time  when  field  energy  reaches  its  first  peak  after  growth. 

In  Fig.  6,  the  total  electrostatic  energy,  ESE,  in  a typical  run  is  shown 
providing  a growth  of  one  to  two  decades  (i.e.,  a factor  of  10  to  100)  in  a 
few  cyclotron  periods. 

(1)  Check  With  Nonlinear  Frequency  Shift  Theory 

One  possible  saturation  mechanism  for  the  drift-cyclotron 
mode  is  the  nonlinear  ion  cyclotron  frequency  shift  [2,9],  which  proposes 
the  scaling  of  the  saturated  electrostatic  potential  as  c,  <*  <65.  We  con- 
ducted a series  of  parameter  studies  with  single-mode  runs  (as  in  the  above 


j 


nonlinear  theory).  The  theoretical  power  dependence,  c “ 6^  is  roughly 

recovered  in  a certain  range  of  6 in  these  simulation  runs  (see  Fig.  7A) , 

where  all  parameters  other  than  6 are  kept  the  same.  The  change  in  6 

comes  either  from  mass  ratio  mj/m^  or  from  relative  plasma  density 
2 2 - 1 

w.j/Wpj,  The  new  scaling  £ « 6 for  6 >0.07  (see  Fig.  7A)  has  no  theore- 
tical origin.  The  predicted  scaling,  c«k,  is  verified  in  a certain  range 

of  tc  in  Fig.  7B,  with  simulation  runs  inwhich  only  k is  varied.  The  other 
-3/4 

scaling,  c * < , for  k>0.7,  has  close  relation  with  a trapping  theory 

and  is  discussed  in  a later  Section. 

In  these  simulation  runs  for  the  parameter  study,  the  real  fre- 
quencies could  not  be  measured  clearly;  thus,  the  classification  of  the 
simulation  mode  as  to  whether  marginally  stable  or  not  as  in  [9]  was  not 
quite  clear. 

Related  to  this  difficulty  of  measuring  real  frequencies  is  the 
generation  of  a zero-frequency  (vortex)  mode  excited  possibly  by  nonlinear 
effects  or  by  quasi  I inear  spatial  diffusion  due  to  charge  separation 
studied  by  Sperling  [10]. 

(2)  Trapping  Theory 


The  phase  space  pictures  in  Fig.  8 reveal  vortex- 1 ike 
structure  after  saturation.  Since  only  mode  3 is  kept,  three  vortices  are 
seen  here  (especially  in  electron  phase  space  pictures).  Also  observed 
in  our  simulations  are  relatively  large  electric  potentials,  which  are 
necessary  (bui  not  sufficient)  for  trapping,  i.e.,  q ^0.3  -0.5  at  the 
saturation  time.  (These  potentials  are  likely  to  be  reduced  with  the  use 
of  larger  mass  ratios.)  This  observation  leads  us  to  examining  the  possi- 
bility of  having  the  familiar  particle  (electrons  or  ions)  trapping  as  the 
saturation  mechanism. 
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We  have  attempted  to  measure  the  fall  (or  trapping)  time, 

Tf,  which  is  the  time  from  the  first  peak  to  the  first  minimum  of  electro- 
static potential,  and  the  potential  at  its  first  peak,  ♦ , in  various 
simulation  runs  to  check  whether  the  parameter  scales  as  as 

expected  when  trapping  is  the  saturation  mechanism  as  in  [11].  This 
attempt  has  not  been  very  fruitful  because  the  measurement  of  and 
carried  somewhat  large  error  bars  due  to  complicated  patterns  shown  In  the 
simulations;  however,  a qualitative  verification  of  the  above  scaling  was 
recovered  from  some  of  our  simulations  (not  shown  here). 

The  discussion  of  electron  trapping  for  a drift-cone  mode  by 

Aamodt  [12]  is  found  to  be  also  applicable  to  the  drift-cyclotron  mode. 

x -vk 

His  approximate  model  provides  the  scaling,  C « <$  x , which  is  his  Eq.  (6). 
Our  simulation  runs  for  6 $0.07  in  Fig.  7A  verifies  the  6^  dependence  (al- 
though this  weak  dependence  could  also  be  as  discussed  in  an  earlier 
“3/k 

Section),  and  the  < 3 dependence  is  more  clearly  recovered  for  k>0.7  in 
Fig.  7B.  Even  without  using  any  approximate  model  (including  the  ion  dis- 
tribution), the  saturation  amplitude  is  by  his  Eq.  (4), 

;.JL  M 1 

“ci  kal  al 
2 2 

where  <>(x)  - ♦Qexp(-x  /A  ) is  assumed  as  in  his  Eq.  (A5A) . This  radial 
structure  is  similar  to  simulation  results  (Fig.  9). 

(3)  Density  Profile  Modification 

As  the  instability  develops  into  the  large  amplitude  stage, 
it  is  observed  as  in  Fig.  10  that  there  is  appreciable  collisionless  (or 
anomalous)  plasma  transport  from  the  high  to  the  low  density  region  due  to 
E * B drifts.  Thus,  the  density  profiles  (integrated  over  the  y-direction) 
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of  electrons  and  Ions  change  their  shapes,  flattening  macroscop i ca 1 ly 
toward  reduced  density  gradients,  which  nay  be  considered  the  saturation 
mechanism  as  discussed  in  Lee  et  al . 1 1 3 J . Again,  this  charge  separation 
in  the  direction  of  density  gradient  might  be  related  to  the  generation  of 
the  col  I i s ion  less  vortex  mode  of  Sperling  {10]. 

in  contrast  to  Fig.  10A,  Fig.  10B  shows  less  modification  of  the 
density  profile  since  this  case  seems  to  have  different  saturation  mech- 
anism as  discussed  in  Sec.  IV  (a). 

(4)  Electron  Heating 

As  the  instability  grows,  electrons  behave  somewhat  coherently 
with  the  wave  especially  before  saturation  (cf,  Fig.  8A)  and  then  develop 
thermal  spreads  near  or  after  the  saturation  time.  Our  examination  shows 
that  this  electron  thermal  spread  is  not  sufficient  to  cause  saturation  of 
the  instability.  The  threshold  electron  temperature  predicted  either  by 
the  local  ROOTS  theory  [3]  or  by  Gary  and  Sanderson  (14)  is  T^/Tj  *6  , while 
the  ratio  at  saturation  in  our  simulation  is  0.07,  thus  eliminating  the 
saturation  possibility  by  the  electron  heating. 
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B.  LOWER-HYBRID  DRIFT  INSTABILITY  SIMULATIONS  USING  ESI  HYBRID  CODE 
Yu-Jiuan  Chen  (Prof.  C.  K.  Birdsall) 

In  the  last  QPR,  Fig.  1 showed  a high  frequency  component 
as  well  as  the  expected  lower-hybrid  drift  mode,  We  tried  using 

different  mass  ratios  (m./me*900  and  1*00)  to  get  rid  of  the  high  frequency 
response,  however,  these  changes  did  not  affect  it.  This  high  frequency 
component  is  still  a puzzle,  but  we  are  working  on  it. 

C.  FIELD  REVERSED  LAYER  SIMULATIONS  IN  Id,  QUAD1 

D.  S.  Harned  (Prof.  C.  K.  Birdsall) 

A nonlinear  quasineutral  Darwin  Code,  QUAD1 , has  been  written 
using  the  framework  of  the  fully  electromagnetic  code  EMI.  The  purpose 
of  this  code  is  to  provide  insight  into  the  problems  of  matching  the  plasma 
fields  at  the  plasma-vacuum  interface.  This  problem  must  be  resolved 
before  moving  to  a two-dimensional  model  to  study  field-reversed  layers 
surrounded  by  a vacuum  region. 

The  model  uses  ion  particles  and  an  infinitely  conducting 
electron  fluid  background.  In  one  dimension  the  field  solver  is  expressed  as 
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This  model  is  essentially  the  same  as  that  suggested  by  other 
groups  working  with  quasineutral  models  |1,2,3).  The  imposition  of  quasi- 
neutrality appears  to  prevent  the  application  of  the  model  in  regions  of 
low  density  (n-*0).  Therefore,  artificial  methods  of  matching  to  the  vacuum 
field  equations  must  be  used  at  some  point  in  the  plasma  where  the  Courant 

condition  ~ < — breaks  down. 

Ax  vA 

At  present  the  code  exists  in  periodic  and  bounded  (conducting 
wall  boundaries)  forms.  It  is  in  its  preliminary  testing  stages.  The  ini- 
tial version  solved  the  field  equations  with  a leap-frog  method.  However, 
to  eliminate  non-physical  high  frequency  effects  this  has  been  changed  to  a 
predictor-corrector  scheme. 
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D.  CONTROL  OF  UNWANTED  BEAMING  INSTABILITIES 
Yu-Jiuan  Chen  (Prof.  C.  K.  Birdsall) 

The  quiet  start  loader  In  ESI  used  to  study  the  lower-hybrid 
drift  instability  has  a physical  multibeam  instability  which  may  affect  the 
simulation  results.  As  the  same  problem  occurs  in  other  simulations  some 
effort  has  been  made  to  study  the  multibeam  instability.  Results  in  Gitomer 
and  Adams  [1]  have  been  reproduced  exactly.  Due  to  the  analogy  between  the 
mult  I beam  and  two- stream  instabilities,  we  tried  to  stabilize  the  multi  beam 
instability  by  adding  a small  thermal  spread  VT1  to  each  beam.  As  shown  in 
Fig.  1 (A,  B and  C) , associated  with  a larger  thermal  spread  is  both  a 
higher  ini  tial  electric  field  energy,  and  a shorter  period  linear  growth 
before  the  field  energy  reaches  saturation.  It  is  noticed  that  all  the 
saturation  levels  are  the  same  as  that  in  the  randomly  Maxwellian  case.  A 
reason  that  a warm  multibeaming  plasma  system  will  reach  the  thermal  equili- 
brium faster  is  that  each  beam  already  is  randomly  Maxwellian.  The  para- 
meters we  used  are  as  follows:  time  step  DT  ■ .04t  , cell  size  x*X_, 

p D 

total  plasma  length  L»128Xp,  number  of  large  groups  NLG-128,  number  of 
particles  n*  16384,  and  the  thermal  spread  added  to  each  beam  VT1»0.0, 

1 OV,  and  0.5  DV,  respectively.  is  the  plasma  oscillation  period. 

The  DV  in  ESI  code  is  defined  as  DV-5.  ,v  VT2  / ( N - 1 ). 
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Section  II I 

CODE  DEVELOPMENT  and  MAINTENANCE 


A.  ESI  CODE 

A.  Friedman 

Bruce  Langdon  of  LLL  has  devloped  an  updated  version  of  the 
ESI  code  which  can  now  be  run  on  either  the  7600  or  the  Cray  using  the  same 
source  file.  This  version,  also  containing  other  miscellaneous  changes, 
is  documented  in  ES1REP,  which  can  be  accessed  by  typing: 


FILEM  RDS  339650  .ESI  ESI L IB  /TV 
ESI L I B ES1REP  OR.  /TV 
TRIX  AC  / T V 

PRINT (<USC  CCC  ESI  REP  BOX  BNN>) 

END 


[trix  print  doesn't  yet 
like  USC  UCB] 


See  also  the  LIBRIS  writeup  for  ESI. 

Stephen  Au-Yeung  and  Alex  Friedman  have  put  together  a ver- 
sion of  the  ZED-compatible  code  ESZ  (for  the  7600  only)  which  reinstates 
the  "original"  history  plots  of  ESI.  The  code  is  called  ESZ1  and  resides 
in  Filem  directory  SAY76  of  user  number  1222.  Changes  to  the  source  are 
flagged  by  "CA"  in  the  first  column  (changes  to  make  ESZ  itself  use  "CADD") . 
Variables  associated  with  the  "original"  plots  typically  have  the  same  names 
as  in  the  original  ESI  with  "0L"  appended.  The  input  variable  "IFOL"  can 
be  set  to  zero  to  suppress  these  plots  if  they  are  not  wanted. 

B.  EMI  CODE 

No  special  report. 

C.  EZOHAR  CODE 


No  special  report. 


D. 


ESI  ♦ EFL  CODE 


[ 

I 


No  special  report. 

E.  RINGHYBRID  COOE 
Alex  Friedman 

This  code  is  under  active  development  and  testing  on  ion  ring 
and  layer  problems.  We  have  had  some  success  In  creating  equilibria  corres- 
ponding to  high-beta  mirror  plasmas,  with  plasma  radius  two  or  three  times 
the  nominal  gyroradius,  but  these  equilibria  exhibit  large  fluctuations  and 
are  thus  too  "noisy"  for  use.  More  flexibility  has  been  built  into  the 
particle  Injector,  and  with  experience  we  expect  to  be  able  to  generate 
reasonably  quiet  equilibria.  We  have  decided  to  apply  the  code  initially 
In  its  present  form  to  field-reversed  mirror  configurations,  in  order  to 
gain  experience  before  attempting  to  modify  the  algorithm  in  the  manner 
described  in  the  previous  QPR.  The  modifications  will  eventually  facilitate 
the  making  of  runs  not  incorporating  the  cold  fluid  ion  component  now 
assumed  to  be  present. 

F.  RJET  DEVELOPMENT 
A.  Friedman 

The  UCB  RJET  is  up  ond  running,  complete  with  the  Versotec  printer 
as  an  independently  controlled  output  device.  Our  initial  impressions 
are  very  favorable,  uith  turnaround  time  and  output  quality  good-  Ue  are 
grateful  to  Harry  Massaro  and  everyone  else  who  helped  get  the  system 
working. 

Some  notes  on  usage  follow. 

To  send  listings,  print  files,  etc.  to  the  Versatec  printer: 

NETOUT  CUCBJ  PRINTF1LELIST  [OPTIONS!  BOX  BNN  ID  / T V 

Options:  8.  8-bit  Ascii  file 

L.  file  contains  lines  longer  than  120  characters 
S.  odd  sequence  numbers  (6-bit  Ascii  files  only) 

ULC.  6-bit  upper /lower  case  file 
(there  are  other  options  as  well) 


i 

\ 


* 
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aUCB*  is  optional  if  one  is  logged  onto  a user  number  uith  UCB  as  the 
default  destination.  Instead  of  NETOUT  one  can  use  BANNER,  which 
has  a wider  range  of  options,  including  *ccsp. * ond  *seq. ■ from  ALLOUT. 

To  send  plot  files  to  the  Versatec  printer: 

NETPLOT  [UCB]  PLOTFILELIST  [OPTIONS]  BOX  BNN  ID  / T V 

Options:  F.  necessary  for  FR80  format  plot  files. 

L.  for  *high-resolut ion*  (I0*xl0*)  plots  - 
default  is  5*x5*.  The  larger  plots  do 
not  fit  entirely  on  one  8.5*xll*  page. 

(there  are  other  options  os  well) 

To  query  the  progress  of  files  through  the  netuorK: 

(CTRL-E)RP?  UCB  and/or  YLU8  UCB  ✓ T V 

To  send  selected  frames  of  a plot  file  to  the  Vdrsatec  while  viewing 
the  frames  on  the  TeKtronix  terminal,  using  utility  routine 
TEKPLOT  (for  DD80  files)  or  FR80PLOT  (for  FR80  files),  type  *USC*. 
and  when  the  routine  asks  for  the  destination  reply  "UCB*. 

To  extract  selected  frames  of  a DD80  plot  Mle  without  having  to  view 
the  file  (and  wait  for  TEKPLOT  to  draw  the  frames),  use  the  utility 
routine  DDEX.  described  elsewhere  in  this  report.  It  is  helpful 
to  have  your  program  generate  a frame  index  to  facilitate  use  of  this 
utility.  Large  plot  files  can  be  sent  to  microfiche,  and  selected  frames 
ser.l  to  the  Versatec:  this  will  ovoid  tying  up  the  printer  for  excessive 
periods  of  time  if  your  plot  file  consists  of  hundreds  of  frames. 


S 


Control  of  the  printer: 

Using  the  operator's  console  (the  teletype  by  the  pdp-11).  one  has 
some  control  over  the  printer  and  the  print  queue.  Uhen  a file  begins 
to  list,  the  teletype  will  print  its  name,  the  word  'begin*,  and  a sequence 
number. 

To  Kill  a printout:  PRNT  KILL  eee.  where  eee  is  the  sequence  number 

(leading  zeros  can  be  omitted) 

To  re-queue  it:  PRNT  REQU  eee 

To  take  the  printer 

offline,  in  order 

to  add  paper,  etc.:  PRNT  DOWN 


To  bring  the  printer 
back  onl i no: 


’ 


ll 
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PRNT  UP 


DDEX  Utility  Routine 

Steve  Au-Yeung  and  Alex  Friedman 


G • 


dd80  Frame  Extraction  Routine 

Availability:  NMFECC  a-machine 

filem  rds  1222  .say 76  ddex 

ddex  allows  the  user  to  extract  selected  frames  from  one  or 
more  dd80  graphics  files  and  write  them  to  a smaller  dd80  output  file  which 
can  be  processed  in  a wide  variety  of  ways.  The  intended  purpose  was  to 
facilitate  use  of  the  "high-quality  hardcopy"  option  of  the  utility  routine 
FROG.  Since  this  form  of  output  is  costly,  large  amounts  of  graphical  out- 
put should  not  normally  be  processed  in  this  manner.  In  addition,  the 
shorter  output  files  (1)  may  be  used  to  create  3S  mm  negatives  suitable 
for  white-on-black  projection  via  give  to  user  - 1 (J  frame)  or  FROG  (full 
frame),  (2)  may  be  sent  to  the  user's  use  versatec  printer  via  netplot, 
or  (c)  may  be  sent  to  the  nips  processor  via  nipit. 

If  your  program  generates  fr80  files  to  obtain  a dd80  file 
use  the  utility  routine  fr80edit.  This  is  especially  useful  if  your  code 
runs  on  the  Cray  which  supports  only  fr80  files.  Send  the  fr80  files  to 
the  7600  via  "netout  a filenames  b." 

Usage:  the  main  usage  form  is  exemplified  by: 

user:  ddex  /tv 

routine/user:  infile:  dxmyplotsl 

routine/user:  box  and  id:  box  bnn  plots1_ex 

routine/user:  > k 17  35  20 

routine/user:  > end 

routine:  all  done 

The  id  string  consists  of  up  to  18  alphanumeric  characters. 
Embedded  spaces  may  appear  incorrectly  and  are  not  recommended.  Use  "_" 
instead.  Note  that  the  frames  need  not  be  in  order  as  they  are  ordered 


- 33  - 


before  extraction.  The  output  file  name  is  constructed  from  the  id  line; 
it  is  of  the  form  dx  [first  eight  characters  of  id).  If  id  consists  of  fewer 
than  eight  characters,  the  filename  is  right-filled  with  asterisks  to 
ensure  a ten-character  file  name.  If  the  file  to  be  created  already  exists, 
the  user  is  given  the  options  of  overwriting  it  or  entering  a new  id  line. 

Some  shorter  usage  forms  are  allowed;  in  particular,  the  input  file  name 
may  be  entered  on  the  execute  line,  or  after  the  command  "infile". 

ddex  is  actually  a "lib"  file;  the  sources  and  binaries  are 
contained  within  the  file  itself,  as  is  this  documentation. 

Options  (all  may  be  abbreviated  by  first  letter): 

thru  extract  all  frames  between  preceding  and  following 
entries,  inclusive. 

infile  change  the  input  file.  All  selected  frames  of  the 
previously  opened  file  are  extracted  before  opening 
the  new  file. 

end  terminate,  processing  all  selected  frames. 

quit  like  end,  but  when  entered  after  infile  command 
causes  previously  entered  frames  to  be  extracted 
without  changing  the  input  file  or  causing  ter- 
mination; more  frame  numbers  can  then  be  entered, 
(useful  when  an  "i"  has  inadvertently  been  entered). 

nochar  eliminates  all  alphanumeric  characters  from  all 

frames.  Useful  for  suppressing  crowded  labels  on 
axes  when  preparing  figures  for  publication. 

offset  set  offset  for  calculation  of  frame  numbers.  Useful, 
for  example,  if  the  user  is  referring  to  a previously 
generated  "nipit"  output  to  determine  frame  numbers; 
tekplot  considers  frame  0 to  be  blank,  and  frames 
1 and  2 to  contain  header  information,  so  that  for 
tekplot  (and  ddex)  plot  information  normally  begins 
in  frame  3.  On  the  nips  output,  frame  1 contains 
plot  information,  so  the  user  may  enter  "offset  2" 
and  then  the  nips  frame  numbers. 

Restrictions:  at  present  ddex  does  not  accept  file  families. 

If  a frame  overlaps  the  end  of  one  file  and  the  beginning  of  the  next,  the 


rout ine  will  fail. 


r 


I 


Support:  this  routine  was  developed,  and  is  supported,  by 
Stephen  Au-Yeung  and  Alex  Friedman  at  the  University  of  California,  Berkeley. 

Inquiries  / gripes  should  be  directed  to  either  of  the  above  at  (A15)  6^2- 1297 
or  6A2-3**77. 

Section  IV 

PLASMA  SIMULATION  TEXT 
Mo  special  report. 

Section  V 

SUMMARY  OF  REPORTS,  TALKS  AND  PUBLICATIONS 
IN  THE  PAST  QUARTER 

No  special  report. 

Eviata. 

In  the  QPR  for  the  previous  quarter: 

Pa9e  $ the  footnote's  author  should  be  Davidson 

page  16  Section  B,  Line  2:  to  read  ..."  x and  k " 

y 

i ' 

H 
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